Disentangling Spatial and Flow Anisotropy 



Boris Tomasik 
The Niels Bohr Institute, Blegdamsvej 17, 
2100 Copenhagen 0, Denmark 

February 17, 2005 



Abstract 

Two generalisations of the blast-wave model to non-central nuclear 
collisions are constructed, and elliptic flow as well as azimuthal depen- 
dence of correlation radii are calculated. Particular attention is paid to 
how different azimuthal dependences of transverse flow direction can cause 
qualitatively different anisotropic fireballs to give same V2 as a function 
of the transverse momentum. The simultaneous dependence of V2 and 
the oscillation of correlation radii on both spatial and flow anisotropy is 
studied in great detail. 

1 Introduction 

The so-called "elliptic flow" V2 [T], observed in non-central nuclear collisions at 
highest available energies has had important implications on understand- 

ing of the collision dynamics in framework of hydrodynamic and cascade models 
^ |2j ■ After exclusion of non-flow effects the "elliptic flow" can be caused by 
azimuthal anisotropy in transverse expansion and/or anisotropic shape of the 
fireball. It was noticed that it was not possible to conclude on the spatial 
anisotropy of the freeze-out state simply from data on V2- However, a con- 
jecture [Sj was made that measurements of HBT correlation radii as functions 
of azimuthal angle |7] will give access to the spatial shape of the fireball. In- 
deed, this was observed |S] in hydrodynamic simulations with two different sets 
of initial conditions: they lead to different final states which could have been 
distinguished by pion interferometry. 

This paper focuses in great detail on the question to what extent the spatial 
shape and the anisotropy of transverse flow can be identified from data. In 
contrast to hydrodynamic simulations where given model and initial conditions 
lead to a single freeze-out state, parameterisations of the final state will be em- 
ployed here. These parameterisations, which will be constructed by generalising 
the blast- wave model [5] , allow to investigate a broad range of various freeze-out 
states and find their possible signatures in the data. A much more systematic 
study than in a hydrodynamic simulation is thus possible; the price to pay is 
that no connection to fireball evolution is made. 
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The generalisation of the blast-wave model to non-central collisions is not 
unique. In order to explore possible ambiguities due to various angular depen- 
dences of the expansion velocity, I will study two models. 

There is a correlation between the spatial and the flow anisotropy in de- 
termining V2', same V2 can be caused by many different combinations of the 
two anisotropics. This correlation depends on the mass of particles but — 
unfortunately — it also depends crucially on the used model. Thus the flow 
anisotropy cannot be disentangled from the spatial anisotropy unless the model 
is known. On the other hand, irrespective the model, the azimuthal dependence 
of correlation radii seems to be mostly sensitive to the spatial anisotropy, at 
least in the low-pt region. 

First, in the next Section I will introduce a generalised blast- wave model. 
Then, V2 (Section O and correlation radii (Sectional are calculated and their 
dependence on spatial and flow anisotropy is studied. Technical details concern- 
ing the calculations can be found in Appendices. 



2 A generalisation of the blast-wave model 

This generalisation follows mainly ref. jlUj with some variations in introducing 
the azimuthal dependence of the transverse velocity. (Experts in the field can 
skip most of this section and look just at the introduction of two different 
azimuthal dependences of the transverse flow velocity, after eq. ((Tn|) .) 

It is assumed that at the end of its evolution the fireball is in a state of local 
thermal equilibrium characterised by a temperature T. Decoupling of particles 
is (almost) instantaneous and can be modelled by the Cooper-Frye formalism 
[TT] along a freeze-out hyper-surface. 

The time of freeze-out does not depend on position in direction transverse to 
the beam, only the longitudinal coordinate matters. Motivated by the Bjorken 
boost-invariant longitudinal expansion jl2| , the freeze-out hyper-surface is given 
by a hyperbola 

To = \? — 7? = const , (1) 

where t and z are temporal and longitudinal coordinate, respectively. We will 
allow for some smearing of the freeze-out time tq by amount Ar. 

We assume that the decoupling matter is distributed uniformly and the 
transverse cross-section has ellipsoidal shape. Thus the density will be propor- 
tional to 8(1 — f) where 

fx2 ^ 

In this equation, x and y are Cartesian coordinates in the direction of the impact 
parameter and perpendicular to the reaction plane, respectively. They can be 
rewritten with the help of the usual radial coordinates 

X = r cos , (3a) 
y = r sin 0s , (3b) 
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Table 1: Summary of model parameters. 



freeze-out temperature T 

average transverse flow gradient po 

variation of the flow gradient p2 

average transverse radius R 

spatial anisotropy a 

mean Bjorken lifetime tq 

freeze-out time dispersion Ar 



(the reason for subscript "s" on the angular coordinate will become clear later). 
The radii Rx and Ry in eq. (|2l stand for the sizes in the corresponding directions. 
They will be expressed via the average radius R and the spatial anisotropy 
parameter a 

Rx = aR, Ry^R/a. (4) 

Thus a fireball elongated out of the reaction plane corresponds to a < 1, while 
a > 1 stands for an in-plane elongated source. 

We do not assume any geometric limitation in the longitudinal direction, 
therefore the fireball is actually infinite in this direction. We can do this as we 
will be only interested in observables at mid-rapidity in collisions at very high 
energy (at RHIC e.g.) where boost invariance is locally established. The actual 
finiteness of the effective source is established dynamically ^]. Those parts 
of the fireball moving too fast forward or backward cannot emit mid-rapidity 
particles. 

The source is modelled by an emission function. This is the Wigner phase 
space density of particle emission 

„, . ,4 mt cosh(2/ -ri) t dr 

Six, p)dx^ ^-^3 dx dy -j== 

xexp(-(L_^)e(l-f)exp(-^). (5) 

Here, we use space-time rapidity 77 and longitudinal proper time t instead of t 
and z 

t = r cosh 77 , (6a) 
z = T smhrj . (6b) 

Momentum p will be parametrised in terms of rapidity y, transverse momentum 
Pt, transverse mass rrit = \/ -f p\ and azimuthal angle 

p^ = (rrit cosh y, pt cos 0, pt sin </>, mt sinh y) . (7) 

The term mt cosh(?/ — 77) in the emission function comes from the flux of particles 
through an infinitesimal piece of the freeze-out hyperbola: p^da^ '11'. In the 
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Model 2 



Figure 1: Two used models of azimuthal variation of the transverse flow velocity. 
The depicted anisotropics correspond to a < 1 and p2 > 0- 



Boltzmann distribution, energy is taken in the rest frame of the emitting piece 
of the fireball 

E*=p''u^ix), (8) 

where is local collective velocity of the fireball. The use of Boltzmann dis- 
tribution is justified as long as the temperature is not too low and the chemical 
potential (for pious) is small; here we put fj, = 0. 

Velocity field u^{x) describes the collective expansion of the fireball f^. In 
longitudinal direction we assume a boost-invariant expansion which is given by 



tanh rj . 



(9) 



The transverse velocity will be parametrised with the help of transverse rapidity 
P 

v± = tanh p . (10) 

Rapidity p will depend on the position in the transverse plane. We will consider 
two models which will differ in the azimuthal variation of the transverse velocity. 

Model 1. In this model transverse expansion velocity is always directed 
perpendicularly to the surface given by f = const ^U] ■ Its angle with respect to 
the reaction plane is thus 

y 



Arctan- 



(11) 



(see Figure nj. Note that 



ta,n d>h — 



tannic = a tan( 



(12) 
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The magnitude of the transverse rapidity also varies with (j)h 

P = fpo{l + P2cos{2(f>b)), (13) 

where po and p2 are tunable parameters. (Note the slight difference to the 
parametrisation of Retiere and Lisa 10 who write p = f(/3o + P2 cos(20fc)).) 
The velocity field u^{x) is then written as 

— (coshp cosh 77, cos^fc sinhp, sin0h sinhp, coshp sinh?^) . (14) 

Later in the calculation it will be convenient to use coordinates f and 0;, in 
the transverse plane instead of x and y. It is shown in Appendix IXI that 

dxdy — Ji[(f)}j) E? f df d(f)b , (15a) 
Ji{(t)b) = (a^cos^ 0h + a^^sin^ 0h)^^ . (15b) 

Model 2. Here it will be assumed that the transverse velocity is always 
directed radially. Then the angle between transverse velocity and the reaction 
plane coincides with 0s (see eq. Transverse rapidity will be given by 

p = fpo(l + p2cos(2(/)s)). (16) 

The difference to Model 1 is in the use of (t>s instead of cf. equation (|13|l . 
The velocity field is similar to eq. H14|l. except for the replacement (f)b <j>s 

u'^ = (coshp cosh?7, cos(j)s sinhp, sin 0s sinhp, coshp sinh?]) . (17) 

In this case, the appropriate coordinates to use are f and 0s- The Jacobian 
is calculated in Appendix IXI 

dxdy = J2{(t>s)R^fdfd(l)s (18a) 
J2(0s) = (a^^cos^^s +a^sin2 0s)"^ (18b) 

All parameters of the models are summarised in Tabled 



3 Elliptic flow 



The elliptic flow coefficient V2 is introduced through the Fourier decomposition 
of the azimuthal dependence of single-particle spectrum PP. At mid-rapidity 
in symmetric collision systems such a decomposition includes only even cosine 
terms 



d^N 



Pt dpt dy d4 
1 d'^N 

— — (l + 2t;2(pt)cos(20) + 

27r ptdptdy 



,). (19) 
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In this formulation, (j) is the angle between the transverse momentum and the 
reaction plane. The coefficient V2 can thus be calculated as 

CPi{pucj>) cos(20)d0 
Jo Pi(Pt,(P)d(j) 

Single-particle spectrum is obtained from the emission function by integrating 
over the space-time 

/ d^xS{x,p). (21) 



Combining eqs. H20() and H21|l we obtain expressions for V2 in our models; see 
Appendix IbI for details of the calculation. For Model 1 (velocity perpendicular 
to the surface) we obtain 



V2 = 



dffjl'' del,, C0S(2</.,) Ji(0,) Ki (!Iii£2!M^) I2 (^£l£i£M^) 



jl df f J^^ dc^k Jl(0b) Ki (iril£££h£(^) I„ (^ p^sin^r-M ^ 

[Model 1] (22) 

while for Model 2 (radially directed transverse velocity) we have 



V2 



f cos(20,) J2(0.) Ki (!2l££SMf:M) (^a^ini^CnM) 



/o'dff/„'"d^,J2(</>.)Ki(^ 



cosh p{r ,(ps) \ T / Pt siiih p(r, 0s ) 
" T 



[Model 2] (23) 



where Ki, Ig, and I2 are modified Bessel functions. Since we integrate over the 
angle, the only difference between the two results is in the use of the Jacobian 
terms Ji or J2. Moreover, the difference between these two terms is only in 
the replacement a — » a~^. As the anisotropy parameter a does not appear 
anywhere else in relations H22I) and (|23|l . any V2 calculated in one model is equal 
to V2 calculated in the other model under transformation a —>■ a^^. Thus we 
have an analytic example of two models which lead to the same V2 , while one is 
elongated in-plane and the other out-of-plane. This clearly demonstrates that 
there is no possibility to distinguish in-plane source from out-of-plane source just 
by measuring V2 ■ 

Physics reason behind the observation that the two models are "inverse" 
to each other can be deduced from Figure ^ The arrows denote expansion 
velocity and their lengths indicate its magnitude. Both situations in that Figure 
correspond to p2 > and a < 1. In case of Model 1, most arrows point rather in 
the reaction plane (which is taken to be horizontal in this Figure), so the major 
boost effect leading to enhancement of the spectrum happens in this direction. 
For Model 2, a larger part of the flow is directed out of the reaction plane and 
the enhancement of spectra due to the boost is in that direction. Thus Model 
1 would lead to positive V2 while Model 2 — in this setup — to negative W2- This 
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is, of course, just a qualitative argument. The equivalence of the two models 
under replacement a ^ is derived analytically. 

We can therefore calculate V2 just for one of the models; results for the 
second one are then obtained trivially. I will choose Model 1. 

Dependence of V2 on the transverse momentum and particle identity in Model 
1 was thoroughly studied in ^U]. In the data [21 E] for small pt, V2 is positive, 
increases with increasing pt and decreases with growing mass of particles. This 
behaviour is reproduced in Model 1 if a < 1 and > 0. (Loosely speaking, 
one of these conditions may be broken, but not "too much" ; see later when the 
results are shown.) In this parameter region, calculation shows that V2 of heavier 
particles can become negative at low pt and start growing and be positive above 
some value of pt |E1 ■ Such a dip to negative V2 is also observed by STAR 
Collaboration for antiprotons in certain centrality bins but the effect may 
not be statistically significant. 

We will be interested in how the spatial and the flow anisotropics are entan- 
gled in determining V2 for various identified particle species. 

In comparing V2 of different species it turns out to be unwise to use the pt- 
averaged V2. This is because the averages are weighted with the single- particle 
spectra which are not alike for different species: those for heavier particles are 
flatter. Hence, averaging V2{pt) for heavy particles can sometimes lead to larger 
resulting values than the same procedure yields with light particles, although 
the value of V2 at any pt is lower for heavy particles. The reason is that flatter 
spectrum for heavy particles gives stronger weight to larger V2 at higher pt. 

Therefore, we study the entanglement of a and p2 in determining V2 for pions 
and protons at two fixed values of pt- As can be seen from eqs. and I^Hjl . V2 
does not depend on i?, tq and Ar. Since the dependence on all other parameters 
was studied in detail in |TO], here we just fix T and po, and plot V2 as a function 
of a and p2 in Figure [5] We clearly see how the two anisotropy parameters are 
correlated and that the correlation depends strongly on the type of particles. 
Recall that a figure for Model 2 would be obtained just by inverting the a-scale. 

Temperature and radial flow can be roughly obtained from azimuthally inte- 
grated spectra which are weakly sensitive to a and p2 . Then, if one was able 
to determine the correct model for the description of the freeze-out, spatial and 
flow anisotropy could be disentangled from measurements of V2 as a function of 
Pt for different identified particle species (Sj. However, the choice of the model 
cannot be based on V2 measurement. 

4 Azimuthal angle dependence of correlation radii 

In this section we focus on azimuthal angle dependence of correlation radii due 
to spatial and flow anisotropy. First, explicit and implicit (/)-dependences are 
discussed. Then, necessary formalism is introduced. Experts can skip this and 
proceed directly to Section lOI where the results are presented. 
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Figure 2: Contour plots of vi as a function of a and p2, calculated in Model 1 
for T = 100 MeV and po = 0.88. Upper row: pions, lower row: protons. Left 
column: pt = 200MeV/c, right column: pt = 500MeV/c. The thickest lines 
show where vi vanishes. Consecutive lines correspond to increments/decrements 
by 0.02. 



8 



4.1 Explicit and implicit azimuthal angle dependence 

In non-central collisions one can study Bose-Einstein correlations of identical 
pions for particles emitted under different azimuthal angles [Tj- To start the 
argument, let us just take a single emitter which emits particles in all directions. 
The directions in which the sizes of the emitter are measured, are given by the 
momentum. Therefore, by changing azimuthal angle (j) of the momentum, the 
correlation radii measure the size of the anisotropic source in different directions. 
This leads to explicit dependence of the correlation radii on i^. (See below 
how correlation radii are defined.) The explicit azimuthal dependence is thus 
connected with the spatial anisotropy. 

In a real case, we have an expanding fireball of which only a part — the 
homogeneity region — effectively produces particles with a given momentum |13| . 
Thus particles in different directions can be produced from different homogeneity 
regions which differ in sizes. This mechanism leads to an additional, so-called 
implicit azimuthal dependence of the correlation radii. It is intimately connected 
with transverse expansion and its anisotropy. Here we want to see how these 
two kinds of effects act together in azimuthally sensitive correlation studies. 

4.2 Formalism 

We will confine ourselves to theoretical calculations at mid-rapidity for sym- 
metric collision systems. The reader is referred to |15l 116) for summary of the 
formalism of Bose-Einstein interferometry in non-central collisions. 

Correlation radii are width parameters of a Gaussian parametrisation of the 
measured correlation function 

C{q,K) = 1 + A eM'Rl{K)ql - Rl{K)ql - Rf{K)qf 

- 2RUK)qoqs - 2Rli{K)q,qi ~ 2Rl{K)q,qi) (24) 

where the momenta of the pair have been parametrised in terms of 

q = Pi~P2 (25a) 
K - i(pi+P2) (25b) 

and the phenomenological parameter A < 1 is due to a variety of effects rang- 
ing from partial coherence of the source up till particle misidentification. The 
standard out-side-long coordinate system is used, with longitudinal axis in beam 
direction, outward axis parallel to the transverse component of K, and sideward 
direction perpendicular to the previous two. Correlation radii are given by sizes 
in these directions. Recall that in non-central collisions we identify the Carte- 
sian x-y-z frame with the collision geometry: z-axis points in beam direction, 
X-axis is parallel to the impact parameter, and y-axis is perpendicular to the 
reaction plane. Hence, there is an angle (j) between the x-axis and the outward 
direction and we are interested in the ^-dependence of the correlation radii. 

If there is no tilt of the fireball in the reaction plane ^T], the two radii i?^; 
and i?^; vanish. This is the case with the used models. The remaining radii can 
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be calculated as [7| 




W^) + (y'» - K(y') - (x2))cos20~ (iy)sin20, (26a) 
W^) + (y'» + \m - (5c2))cos20+ (iy) sin 20 



'OS 



+ - 2A(tx) cos (/) - 2/3t(ty) sine/) , (26b) 

i((y2) - (x2^)sin20+ (xy) cos 20 





where 



{f{^)){K) 



U{x)S{x,K)d\ 

^ S{x,K)d'^x 
x^ - (x^) . 



X 



(27b) 



(27a) 



The explicit azimuthal dependence is displayed in eqs. (|26l) . In addition, the (co- 
)variances (x^x'') can depend on (p and this is the implicit azimuthal dependence. 
From eqs. (|26|) . correlation radii can be calculated for both Models just by 
inserting the corresponding emission function. 

Azimuthal dependence of the correlation radii can be analysed with the help 
of Fourier decomposition. Due to a number of symmetry arguments ^,16, in our 
setup, the relevant Fourier series, truncated after the leading oscillating terms, 
are 



4.3 Results 

We will focus on and i?^, as we are interested in anisotropics in the transverse 
plane. The absolute sizes of these radii together with their oscillation amplitudes 
scale with the total geometric size R. We can get rid of this scaling and thus 
observe the effect due to anisotropics more cleanly when wc study the ratios 



These ratios for Model 1 and Model 2 are plotted in Figure|2| In most cases, 
oscillations of correlation radii are mainly determined by the spatial anisotropy 
and not so much by the flow anisotropy. The only exception is i?^ at high pt 
in Model 2: the 0-dependence in this model changes from shape-determined to 
flow-determined, i.e. dominated by the implicit azimuthal dependence. In all 
other cases, 0-dependence of the correlation radii follows the explicit azimuthal 
angle dependence rather well. 

The behaviour of i?^ in Model 2 is shown in Figure 01 The azimuthal angle 
dependence changes qualitatively with increasing pt- at the given parameter 



Rli^) 
R'si'^) 
Rlsi^) 
RfW 



<o + 2<2 cos 20+... 
i?f + 2i?2 2cos20+ . . . 
2i?L2sin20+... 



Rf o + 2i?f 2 cos 20 ... . 



(28a) 
(28b) 
(28c) 
(28d) 



Rh/Rlo and RI2/RI0 M- 
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values oscillation amplitude is positive at low pt and becomes negative at large 
Pt- 

Now we come to the question if one of the models can be disqualified by 
comparing to data. In Figure [51 we see two models which reproduce V2 quite 
well (this is shown in JQ); they are related by transformation a ^ a~^. How- 
ever, since the oscillation of correlation radii is mostly shaped by the spatial 
anisotropy, the two models lead to opposite predictions of the sign of the oscil- 
lation amplitudes — and Model 2 is in qualitative disagreement with the data. I 
do not try to find the perfect fit here; this turns out to be a problematic task 
with the blast-wave model |2U[ I18| . Nonetheless, note that the qualitative fea- 
tures of the azimuthal angle dependence of the correlation radii are reproduced 
in Model 1 under the assumption of an out-of-plane extended source, which 
confirms the earlier conclusion of STAR 

5 Conclusions 

Generalisation of the popular blast-wave model to non-central collisions is not 
unique. Many possible ways differ in how the transverse velocity depends on 
the azimuthal angle (ps- 

In this paper two such generalisations were constructed. Then, a number 
of statements scattered in literature were demonstrated in a unified framework 
of a generalised blast-wave model, which is often used in many variations. The 
interplay of HBT analysis with V2 for identified species should be stressed. 

I showed analytically how two very different fireballs can lead to the same 
V2, such that from measuring only this quantity one cannot conclude whether 
the source is elongated in-plane or out-of-plane [Hj. 

The azimuthal angle dependence of the correlation radii, on the other hand, 
can be used for this. It is mostly sensitive at low pt to the spatial anisotropy of 
the fireball 0. 

When the type of the model is identified from comparison to data on two-pion 
correlations, spatial and flow anisotropy can be disentangled from reproducing 
V2iPt) of different identified species and (j)- dependence of correlation radii. Before 
measuring anisotropics, temperature and radial flow can be determined from 
azimuthally integrated single-particle spectra which are nearly independent of 
the anisotropy parameters. 

Among the two models used here. Model 2 fails qualitatively in reproducing 
data on azimuthal angle dependence of the correlation radii. In a simple qual- 
itative comparison with the data, the other Model indicates that the observed 
fireball in non-central Au-|-Au collisions at -^/s = 200 AGeV is elongated out of 
the reaction plane, in agreement with conclusions of jl9| . 
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A Jacobians for integration in transverse plane 

As it is simpler, we begin with the Model 2. The aim is to use f and (ps as 
coordinates in the transverse plane. Because 

X = r cos 0s (29a) 

y = r sin 0s, (29b) 

we only have to replace r by f. From eqs. ^ and Q| we obtain 

r 



This leads to 



, a^sin^^s. (30) 



f R cos ( 



\J cos^ 0s + o? sin^ ( 

f R sin 0s 
\J cos^ 0s + a? SIT? i 



(31a) 
(31b) 



and 



, , R'^rdrd(j)s , , 

= ~2 2J. , 2 ■ 2 . • (32) 

a cos^ 0s + sm 0s 



Thus we derived eq. ((TS|l . 

For the Model 1 we want to use the angle 0h as a coordinate instead of i 
By making use of eq. H12|l we can rewrite eq. (|31|) into 



fi? sgn(cos 0s) 
i/a^2 + tan^ 0s 
f Ra^ cos 0b 



y = 



V cos^ 0b -f- 


- a 2 sin^ 0b 


f i? sgn(sin 




a^2(;.otan^0 


■s +a2 


rRa-'^ 


sin 0b 


\/ a? cos^ 0b -f 


- sin^ 0b 



(33a) 



(33b) 



In the last expression we introduced 

, , f +1 : a; > 

'S^(^) = l -1 : x<0 

and exploited that 

sgn(sin0s) = sgn(sin0b 



sgn(cos0s) = sgn(cos0b) 
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From eqs. (|33f) it is straightforward to obtain the Jacobian for Model 1 

d^dy- , f/f't- (34) 

cos^ cpb + a ^ sm (pi, 

Notice that apart from the use of cpb instead of cps , one obtains the Jacobian for 
Model 2 from that of Model 1 just by replacing a a^^. 



B Calculation of V2 

Let us calculate V2 for Model 1. First, we need the azimuthally integrated 
single-particle spectrum in the denominator of eq. H20I) . From eqs. ^ and H14|l 
we obtain that 

p^Wp = mt cosh(?7 - y) coshp(f, 4>b) - Pt sinhp(f, 0b) cos((/) - 0b) . (35) 

This is the energy argument for the Boltzmann distribution. In accord with 
eqs. and ©, azimuthally integrated spectrum is obtained as 



r2Tr 

PliPt,' 



2tt 

4 



d(j) J d X S{x,p) 

rdcPbMM I'drrR' r 



(27r)3 



dT T ( {'T ^ '''a 



exp 



VQjFX^ \ 2 At 



\2 



( mt cosh{r] - y) cosh p{f,4>b) - Pt smhp{r, (j)b) cos{(j) - (j)b)\ 
X exp I 1 (36) 

where Ji{(j)b) was defined in eq. pSII . The integration over r is trivial. We are 
interested in mid-rapidity particles in the centre-of-mass frame, so y = 0. Then, 
integration over rj can be performed and leads to the modified Bessel function 
Ki We can exchange the order of integrations in (p and (pb, and perform 

a transformation — > (p — (pb = ip- The integral in ip can then be performed 
analytically and leads to the modified Bessel function Iq. We finally arrive at 

n Turrit 



Pi{pt,4>)dcp = dffj^ d<PbJi{<Pb) (37) 

I'mt cosh p{f,(pb)\ ^ f Pt smhp{r,(pb) 
xKi — io 



T J \ T 

The numerator of eq. H22|l is obtained in a similar way as the azimuthally 
integrated spectrum, we just add a factor cos(20). After performing the inte- 
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gration over r and rj we obtain 



27r 







Pi{pt,(t)) cos{2(t>) d(f> 

X fj dct> cos(20) exp eos(0 - 0,)) . (38) 

Now again, we write (j) — ip + (l)b and decompose 

cos(20) = cos(2V' + 20f,) 

= cos(2?/') cos(2(/()f,) — sm{2ip) sin(20b) . 

The ^/i-integral with the term proportional to sin(2^/') vanishes. The second 
term, proportional to cos (2-0) exp cos?/') leads to a modified Bessel function 
I2 PH- As a result we thus obtain 

Pi{pt,(t>) cos{2(j>) d(j) 
27r2 











drr d(t)b Jii4>b) cos{2(j)b) 
Jo 

/' m f cosh n(f, (hf,)\ /n* sinh p(f, (?!)!,) \ 

X Ki f ^ Y^^j h i- j . (39) 

By dividing this equation with the denominator derived in eq. H37|l we obtain 
the expression (|22|l for V2. 

Calculation for Model 2 follows exactly the same steps as we've gone with 
Model 1. The only difference is that 4>b is replaced by (f>s and one uses J2(0s) 
instead of Ji(0(,). 
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Fi giire 3: Dependence of the normalised second order oscillation terms i?o^2/ -^o.o 
(upper rows) and P?siI^sQ (lower rows) on spatial anisotropy a and flow 
anisotropy p2, calculated for Kt = 300MeV/c (left columns) and Kt = 
900MeV/c (right columns), with Model 1 (upper panel) and Model 2 (lower 
panel). Thick contour lines correspond to 0, consecutive curves to incre- 
ments/decrements by 0.1. Other model parameters in the calculation were 
T = 0.1 GeV, po = 0.88, R = 9.41 fm, tq = 9 fm/c, and Ar = 1 fm/c. 
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Figure 4: Azimuthal angle dependence of i?^ in Model 2 at various transverse 
momenta. Values of parameters used in the calculation: T = 0.1 GeV, po = 0.88, 
P2 = 0.2, ii = 9.41 £m, a = 0.95, to = 9fm/c, and At = 1 £m/c. Different curves 
correspond from top to bottom to transverse momenta of 0.2, 0.4, 0.6, 0.8, and 
1 GeV/c. 
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Figure 5: Azimuthal angle dependence of correlation radii in Model 1 (left) and 
Model 2 (right) in comparison with data. Curves and data points correspond 
from top to bottom to transverse momenta of 0.2, 0.3, 0.4, and 0.52 GeV/c. 
Values of parameters used in the calculation are chosen such that V2{pt) is 
reproduced (not shown here, see [TH|): T = 0.12 GeV, po 0.99, p2 = 0.035, 
R = 9.41 fm, To = 5.02 fm/c, and Ar = 2.90 fm/c. Spatial anisotropies are 
a = 0.946 (Model 1) or a = 1.057 (Model 2). Data points are measured by 
the STAR collaboration in 20-30% centrality events of Au+Au collisions at 
Vi= 200ylGeV. 
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